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The  ever  increasing  need  to  be  informed  has  caused  an 
ever  increasing  volume  of  publications.  The  very  quantity 
of  these  publications,  however,  makes  it  difficult  to  satis¬ 
fy  the  need,  because  it  becomes  more  and  more  arduous 
to  locate  the  literature  pertinent  to  any  one  subject;  and 
because  any  one  of  us  becomes  more  and  more  weary  of 
having  so  much  to  read.  Moreover,  only  a  small  minority 
of  these  publications  may  help  satisfy  the  need,  because 
an  overwhelming  majority  is  void  of  content. 

Consequently,  it  is  imperative  to  present  to  the  public 
only  manuscripts  which  do  carry  new  information  without 
undue  verbiage;  hoping  to  be  judged  not  by  the  number 
of  published  pages  but  by  their  quality.  In  particular, 
papers  in  the  field  of  Applied  Mathematics  ought  to  be 
published  only  if  they  contain  one  or  more  of  the  follow¬ 
ing  items:  new  basic  results,  new  methods,  new  applicat¬ 
ions,  new  numerical  results,  new  presentation  of  difficult 
and  important  topics,  up-to-date  bibliographies;  and  if 
the  number  of  their  pages  is  not  dictated  by  the  desire  of 
imposing  upon  the  superficial  reader. 

To  discharge  our  contractual  obligations,  we  publish 
Technical  or  Scientific  Reports,  such  as  the  one  you  now 
have  in  your  hands.  It  has  been  our  constant  policy  to 
see  to  it,  that  they  satisfy  the  above  strict  criterion. 
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Electron  Straggling 
Theory  and  Calculation 

Nathan  Grier  Parke  III 


Introduction. 

In  the  course  experimental  work  on  the  effect  of  the  passage 
electrons  through  a  foil  on  its  resistivity^  this  laboratory  undertook 
some  routine  calculations  of  the  distribution  of  the  electron  energy  losses 
about  the  mean  at  different  depths  in  the  foil.  The  initial  reference  used 
was  ”The  Passage  of  Fast  Electrons  Through  Matter;  by  R.  D.  Birkhoff,  in 
Volume  XXXIV  of  the  Encyclopedia  of  Physics,  Springer,  Berlin  1936 
pp.  53“138-  This  starting  point  quickly  lead  to  the  classic  paper  of 
L.  Landau:  J.  Phys.  USSR,  8  201  (19^^)  and  a  paper  by  T.  Blunck  and 
S.  Leisegang:  Z-  Physik  128  500  (1950).  We  had  some  difficulties  making 
a  calculation  that  had  sufficient  validity  and  avoided  enough  complicatio'^ 
to  be  useful  and  tractable.  A  report  of  the  technical  details  of  this 
work  seeas  Justified. 
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Iiaportance  of  Landau's  Theory, 


A  diacuaaion  of  Landau's  paper  and  theory  is  inportant  because  of 
the  role  it  plays  in  the  subsequent  Blunck  and  Leisegang  work  and  in 
the  Birkhoff  article  in  the  current  edition  of  the  massive  Handbuch  of 
Physics^  edited  by  S.  Flugge  and  published  by  Springer^  Berlin.  In 
addition^  there  is  an  article  by  Wolfgang  Borsch-Supan:  "On  the  Evaluation 
of  the  Function 


(1) 
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for  Real  Values  of  ^  Journal  of  Research  of  the  National  Bureau 
of  Standards  -  B.  Mathematics  and  Mathematical  Physics^  65B  2k^-2^0  (I96I). 

This  function  is  Landau's  "universal  function"^  describing  the 

transformation  of  the  spectrical  distribution  of  energy.  It  is  the  function 
that  was  first  calculated  numerically  by  the  calculation  Bureau  of  the 
Mathematical  Institute  of  the  Acadeuqr  of  Sciences  of  the  USSR.  The 
Borsch-Supan  paper  presents  a  table  of  the  function  to  five  signi¬ 

ficant  figures  covering  the  >  range:  -k  (0.CE5)  4.8  (O.l)  20  (0.2)50(0.5)l0 
Blunck  and  Leisegang  base  their  work  on  a  +  2^  approximation  of  the 
Russian  calculations  as  the  sum  of  four  Gaussian  distributions^  i.e.^ 


(2)  ♦ix)  = 


with  coefficients: 
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2 

3 
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0.174 

0.058 

0.019 

0.007 

_ 

0.0 

3.0 

6.5 

11.0 

1.8 

2.0 

3.0 

5.0 

valid  ov«r  tha  range  X  ^  ^  ^  15  ,  with  an  error  of  +  2^. 
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By  "straggling'*  is  meant  the  variation  in  the  amount  of  energy  lost 
by  successive  par  tides  pas  sing  through  the  medium.  To  calculate  the 
straggling^  is,  in  theoretical  terms,  to  find  the  probability  density 
of  energy  loss  ^  after  the  electron  has  traveled  a  distance  X  in  the 
medium.  Landau's  solution  is  related  to  his  universal 

function  in  the  following  way 


wv*  ZA 


where. 


4  ~ 


\  -  >■ 


_ 

C  ^  0.577  ‘  Euler's  constant 


Ui'  •  -U  I  ^  ^  t: 

I  ’1.2 


I  =  II.  i  eV 
0 


an  ionization  potential 


A  Discussion  of  the  Landau  Theory. 


Landau  bases  his  solution  on  differential-integral  equation. 
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vhere  ^  W  is  characterized  as  the  probability  (per  unit 

length  of  path)  of  an  energy  loss  t  for  a  particle  of  energy  where 

the  ionization  losses  are  considered  to  be  soiall  coBipared  with  .  For 
Vsrtf)  he  takes, 


(5) 


W(t)  - 


Tc  M  j  Z  Z  I 


where. 
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Avogadro's  number 
mass  of  electron 
charge  of  electron 
related  to  by 
density  of  substance 

sum  of  the  atomic  numbers  of  the  molecules 
in  the  substance 

sum  of  the  atomic  weights. 


The  reasons  given  for  equations  (4)  and  (5)  left  the  writer  with  an  uneasy 
feeling  and  an  attempt  was  made  to  find  a  more  satisfactory  derivation.  One 
specific  cause  of  uneasiness  is  the  fact  that  the  classical  probability 
density  for  energy  loss  £  "per  collision"  is  ec  l/c^in  fact 


(6) 


W|(i) 


[ 


O  otherwise 


If  there  are  collisions  per  unit  distance,  and 

(7)  i'--  ‘  *t„ 

is  the  random  sum  of  the  energy  losses  in  these  collisslons^  the  t 
is  the  energy  loss  per  unit  distance  and 
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iri  trie  dLStributi.on  of  the  insLe^d  of  (>),  The  only  conclusion  one 
cm  ^-^jach;  knowing  Landau^s  stature  and  ability^  is  that  his  results  are 
correct  c-.*  the  order  of  accuracy  he  has  in  mind  but  their  justification 
^‘ar  be  s^.aced  less  elliptically.  This  can  be  done  if  we  derive  equation  (4) 
“cc‘iii.sion  by  collision'^  rather  than  "per  unit  distance". 

To  begin  with  let  us  compute  the  distribution  function 


f  (1+',  a) 


probability  density  for  energy  loss  ^ 
after  1  collisions 


a  function  of 

probability  density  for  energy  loss  A 
after  n  collisions 


It  IS  clear  that 
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fi,  A 
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<30 
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and  hence  that  the  finite  difference 


oo 


\M 


If  we  write 

Wi  »  f  I  C 


average  collisions  per  unit  length 
h*  hoX  collisions  in  length  X 
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then 

Qo*) 

since. 

o 


(10) 


(y^,  /i  ) 


*1o 


Ho  f  :fcn-vi,ik)  -  f6  ji?.'') 


and  dropping  the  bar  and  writing 

Wet)  =  no  vw,cf) 

oo 

=.  Jwet)  i^-i)  - 
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(10”)  (>/&) 


we  obtain  the  Landau  differential-integral  equation  (5)  and  need  an 
estimate  of  f\^  •  In  the  first  place 


where  b  is  the  "impact  parameter"  of  the  collision.  The  electrons  thac 
will  "collide"  with  the  incoming  electron  will  lie  within  the  range 
and  b^  corresponding  to  max  and  min  energy  loss  per  collision.  For 
an  electron  traveling  unit  distance  the  colliding  electrons  will  fill 
volume  TtCb:-w:yi 


(13) 


C  CiH?  / 


The  density  of  "fixed"  electrons  is 


(14) 


Nfl^ 

ZA 


&lecfrwr>S  \ 
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80  that 

(15) 


•£.  ZA  ( 


CcletW«x^$  \ 

j 


(16)  WCt)»  p. 


where, 

(17) 


ire^  Z  £ 

E«  XA 


The  physical  argument  involves  two  reasonable  asstamptlons ; 

A.  it  makes  no  physical  sense  to  estimate  closer 

than  a  difference  quotient  involving  a  single  collision. 

B.  is  a  sufficiently  good  estimate  of  the  depth  reached 
in  A  collisions  so  that  we  can  equate 

(18)  ^  Cf(''  +  tyA)  - 

The  preceding  discussion  shows  that  £  in  WCl)  stands  strictly  for  the 
energy  loss  per  impact  and  «  UZ/i)  is  the  probability  density  for 

this  loss. 

Landau *8  Approximating  Solution  of  His  Differential  Integral  Equation. 

Landau's  approximating  solution  of  his  differential-integral  equation 
was  reached  by  the  use  of  Laplace  transforms.  Let 

(19)  ^ 
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then,  by  the  Inversion  theorem, 

(20) 


Taking  the  Laplace  transform  of  equation  (5)  gives 


(21) 


because. 


o 

00 

1 


f(R,  li'JJA' 


Equation  (21)  la  now  an  ordinary  differential  equation  in  w  aatlsfylng  the 
initial  condition, 

•|-(o;h.)s  S(a)  or 
The  solution  Is  therefore 

(22)  =  e 
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or,  more  clearly 

(l-  C, 

(22')  ■*  e 


moment  generating  function 
for  the  probability  density 
WtCC)'  for  energy  loss  per 
collision* 


where y 
(23) 


This  Landau's  general  expression  for  the  distribution  function  f  in  terms 
of  the  function  WCt)  where  wCl)  is  the  number  of  collisions  per 
unit  path  length  with  energy  loss  between  1  and  To  proceed  further, 

Landau,  in  effect,  takes: 


(16) 
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To  eee  this,  we  follow  his  steps  where  he  approximates  the  Integral 


(25) 


00  ^  ^ 


where  £i  4  1  ^  •  To  evaluate  the  first  integral  on  the  right 

hand  side  he  refers  to  a  result  from  Livingston  and  Bethe^  Rev.  Mod. 

Fhys.  £245  (1937),  i.e.. 


(26) 


^  twu)at  *  a  ^  a 


We  notice  that  if  we  take  *  ‘T*.  In  the  region  <  1^1 


=»  ft  ^  ~ 


and  that  in  effect  Landau  uses  the  law 

wai  =  Vi‘ 

over  the  interval  to  «  In  addition  he  makes  the  approximation 

in  the  region^  f  1  <  £  <  •  This  will  work  for  sufficiently  small 

Integrating  the  second  term  on  the  right  by  parts 

f  ft- 

Substituting  I  ^  noting  that  again  that 
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oO 


gives 


^  f . 1, 

t  I 


where 


U  + 


otl 

r  -i 

e  A  .  _  e 

2 

I 


1 


I 


(Ryshlk-Grads  tein : ”Tab les 
of  Series^  Products^  and 
Integrals^'*  VEB  Deutscher 
verlag  der  Wlss^  Berlin^ 

1957  p.306) 


where 


Hence 


(27) 


and 


(28) 


C>  >  0.^77216  .,«,Euler'«  constant 


^ ,  t(i  -  C 


f 


WU) 


We  can  now  write. 


f  -  -L  \  e.  ^  «>» 

•*  lire,  J 


If  we  write 

(29) 


f 


^  sax 


4-  - V. 


I 


u=^  f  ^ 

+  d-i 
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Then 

(30) 


where 

(31) 


I 

Z-ri 


-r«o  •♦ff' 


is  Landau's  "universal  function,  carefully  studied  and  calculated  by 
Borsch -Supan. 

Physically  it  is  seen  that  Landau  makes  two  assumptions 


A.  The  distribution  w  £<.)  -m. 


where 


'  4i:'® 


B.  We  don't  have  to  estimate  because  the  contribution  to  the 

energy  loss  is  negligibly  increased  by  taking  f ^  »  *0 

instead  of  f • 

An  Alternative  Approach  *  The  Method  of  Moments. 

After  critically  examining  Landau's  solution^  several  alternatives 
to  obtaining  a  tractable  solution  suggest  themselves.  In  the  first  place^ 
if  W|CC)  is  the  probability  density  for  the  energy  loss  in  a  single 
collision  and  we  take  to  be  the  number  of  collisions  in  going  a  distance 

nien  the  moisent  generating  function  ^  )  is  given  by 

(32) 
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where ^ 

(33) 


oo 


but 


Or  we  get 


(33') 


in  ( I  -  i  ) 


*  *  4 


e 


-h.ic 


m 


to  compare  with  the  Landau  solution  via  the  differential-integral  equation 

4)(x,  j.)  =  e  ‘  -» 

To  compare  the  two  in  another  vay^  e.g.  in  their  estiute  of  the  first 
two  moments^  notice^  that 


C,  (j.1  »  «- 

# 

-  [  w.u,  z 

•  1 


where 


or 


C.Cv)  ’  ^  T-  v‘”  “ 


hence  the  name  moment  generating  function. 
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If 


then 


(r( 


1  - /V t 1: 


^  -A 


=  B  ^  iK, 


^x-t 


^  r>.<  [  <-  K|»  +  (“f<, 


>  If 


V4> 


=  A  ^  -  -n.xi 


*d*4tx,  |») 


r  to  1 

*  hox  •+ 


p-/\^  ^ 
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li  ^  't4>€  otVver 

4c.,p)=  ^  RfV- 

=  ^ i  Cp.p- 


.'50  chat 

/V  =  -  ia,VkC 

B  H-’’  ■*■ 

\  B-  A*"  -  fix 


The  implication  of  this  result  is  that  the  variance  of  ^  computed  by 
the  Landau  differential-integral  equation  will  be  too  large  by  a  factor 
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Hence ^ 


I 


where 


hence ^ 


(-  ..  - - 


■*"  ■»  -i  4iJi®5  ijuo^  •■=.  In  a  case  of  Interest, 

'<■%  T  iu 


i 


I 


1  - 


1. 


A  %>U' 


(  ^  -  t'\  •>  it  *  0-^  =•  4-. 


Also, 


lu^ 


6l-> 


/J, 


&0 


It  would  seem^  as  might  be  expected  on  the  grounds  of  the  reasonableness 
of  both  models ^that  the  numerical  error  is  small.  Nevertheless^  the 
expression  based  on  random  walk^  i.e.^ 


(32) 

where 

(33) 

where 


‘♦>CK,|r)«  [c,c|.^3 


oo 


^i<r)  |^w,a)e“^‘4£ 


WiCC)  distribution  of  f  ,  energy  less  per  collls 


lor. 
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seems  coaceptuaily  simpler  and  may  lead  to  a  more  tractable  inversion  to 
CD tain  1\)  •  This  inversion  is  being  worked  on  but  has  so  far  resisted 

all  attempts. 

On  the  satisfactory  solution  of  the  inversion  problem  rests  the  easy 
calculation  of  straggling  due  to  windows  in  series. 
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